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The improved staggered quark action has enabled breakthrough calculations of QCD with three light dynamical 
quarks. The largest remaining scaling violations come from "taste-changing" interactions, which are signifi- 
cantly reduced over naive quarks because of carefully chosen smearing of the gauge links. Here we examine further 
improvements to the quark action and show that additional smearing of the links can reduce the taste-changing 
errors by another factor of two, whilst retaining a fully a^-improved action. The taste-changing errors can be 
understood within a perturbative calculation of the four-quark operators that appear at one-loop. We also discuss 
unquenching techniques with such a highly improved action. 



Accm'ate lattice simulations of the Standard 
Model require numerical simulations in lattice 
QCD which include three light flavours of dy- 
namical quarks. The improved staggered-quark 
(Asqtad) formalism is the only one capable of 
delivering large numbers of configurations with 
three small quark masses anytime in the near fu- 
ture Recent calculations |21 have shown that 
such a formalism is already highly successful, de- 
livering numbers for "gold-plated" quantities that 
agree with experiment at the few percent level. 
This analysis, while impressive, highlighted that 
increased efficiency in the simulations would be 
possible if the action could be further improved. 
Here we systematically explore this possibility. 

Staggered quarks have 16 non-degenerate pi- 
ous made from the four "tastes" — the remnant 
of the doublers — whose masses do not vanish 
for zero quark mass. The splitting comes from 
mixing between the tastes and vanishes like a?. 
Nai've staggered quarks suffer from poorly con- 
vergent perturbative expressions and large pion 
splittings which are suppressed by a factor of 2-3 
for typical simulations by the use of fat links 
consistent with the reduction of a pure a? ef- 
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feet to agO^. Staggered quarks which are im- 
proved to (^(asa^, a^)-errors 0], significantly re- 
duce splittings in the pion spectrum |2| , and have 
small renormalisations [Sj . A related scheme with 
smaller pion splittings is the HYP action |7j. 

We can further reduce these undesirable split- 
tings by adding explicit 4-quark counter-terms to 
the lattice action [H]. These are computed in 1- 
loop perturbation theory (since taste-changing is 
perturbative) . As we show, a simple alternative is 
to suppress the 4-quark counter-terms by modify- 
ing the quark-gluon vertices. Using insight drawn 
from perturbation theory, we have designed a new 
quark action that has comparable taste-changing 
interactions to HYP, but is much simpler to im- 
plement and unquench and retains the essential 
a^-improvement necessary for accurate simula- 
tions. 

Taste-changing is most easily understood in 
terms of naive quarks, which are equivalent to 
4 identical staggered quarks. Nai've quarks have 
an exact symmetry: 

quark(p ^ 0) = quark(p C'/a) (1) 

where C = (1,0,0,0), (1,1,0,0), ... giving 16 
identical copies or "tastes" of every quark, one in 
each corner of the Brillouin zone. These differ- 
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Figure 1. Sample taste-changing diagrams for 
massless naive quarks 



ent tastes can mix by exchanging a hard gluon 
as shown in figure [r(a)| This gluon is highly vir- 
tual with momentum 0{n/a) and thus the quark- 
quark interaction is effectively a purely pertur- 
bative contact interaction at typical lattice spac 
ings 



The tree-level interaction of figure 1(a) 
was understood and completely removed with 
the improved staggered action by using smear- 
ing to suppress high-momentum gluon emission 
from quarks |3l4j . At one- loop order, there are 
5 non-zero diagrams that contribute to the taste 



changes; one is shown in figure 1(b) In the mass- 
less limit there are 21 independent counter-terms 
coming from 1-loop diagrams: 
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It is not obvious which of these terms is most 
important, but a crude estimate of the taste- 
changing is given by kcl- We want to min- 
imise this quantity when designing new, further 
improved, quark actions. 

One-loop taste-changing requires large momen- 
tum transfers through the gluons, consequently 
we can suppress these contributions by soften- 
ing the quark-giuon vertices at high momentum. 
This softening can be achieved by carefully smear- 
ing all of the link operators in the quark action. 

Extensive smearing leads to a large number of 
terms, many of which involve very long products 
of link operators. While 1-gluon vertices are well 
behaved by design, 2-gluon vertices (tpA^ip) tend 
to become huge and result in large taste-changing 
contributions. This problem can be minimised by 
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Figure 2. A schematic of Fat7 smearing, and 
the principle behind smeared-smearings: the orig- 
inal smeared links are used in the second smear- 
ing. The is because appropriate coefficients and 
sums over all directions of the path elements are 
not shown. 



re-unitarising the smeared- link operators. 

Tree-level taste-changing was removed by 
"Fat 7" smearing (figure I^J 0], but this leaves 
the side links unsmeared. A better choice is to 
apply Fat7 twice: Fat 7(g) Fat7 (figure 2b). To 
avoid the combinatorial growth of ipA^ip ver- 
tices, the links are re-unitarised after each smear- 
ing: Fat7_R(8)Fat7fl,. The a^-improved cousin, 
Fat7/{(g)Asqtad, uses the same re-unitarised Fat7- 
smeared links in the original Asqtad action, 
but with only one re-unitarisation step. We 
also investigated SU(3) back-projection instead 
of re-unitarisation, however, it is much harder 
to unquench for no discernible benefit. We 
have developed two methods to compute the un- 
quenched molecular dynamics force term [Hj. Re- 
unitarisation, Vr = {1/VVV^)V, can be differ- 
entiated via the chain rule ^U] . This can be done 
explicitly or approximated stochastically, the two 
methods give consistent results but the stochastic 
estimation, while introducing additional noise, is 
faster. The Fat7fl(g)Fat7ii; action, like the HYP 
action, has tree-level errors. These are absent 
in the Fat7fl(8) Asqtad, and this action has ~ 2 
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Figure 3. Quenched, /? = 5.93, m = 0.03 simulations with Wilson glue of the splittings in the squared 
masses of 1-, 2- and 3-link pious above the Goldstone. Errors are statistical, and come from uneven 
numbers of configurations. The perturbative estimate of the taste-changing shown on a different scale 
qualitatively tracks the simulation data. The proposed action Fat7i?(X)Asqtad, shown at the right, is 
about two times better than Asqtad, shown on the left. 



times smaller taste-changing than the traditional 
Asqtad action. This suggests that FatTi?® Asqtad 
is the most highly improved action to date. We 
compare the trend in X] kcli computed perturba- 
tively, with quenched simulation results for the 
splitting in squared pion mass in figure |21 The 
perturbative results track the simulation results, 
confirming that taste-changing is essentially per- 
turbative in character. 

These studies confirm that FatTijCgiAsqtad is 
the most accurate light-quark action available, 
with residual discretisation errors of 0{a^as), 
and taste-changing that is substantially smaller 
than for Asqtad and comparable to HYP. 
Fat7fl®Asqtad is far simpler to code than the 
HYP action, more accurate and easier to un- 
quench. Finally, these studies demonstrate that 
taste-changing is essentially a perturbative phe- 
nomenon. Significant further improvements, 
through more complex smearings, are probably 
difficult or impossible. 
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